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ABSTRACT 

Aims. We present ground-based observations of the disk around the A-type star J3 Pictoris to obtain scattered light images at the 
highest angular resolution (60 mas, equivalent to about 1 AU at the distance of the star) and the highest contrast in the very close 
environment of the star The purpose of this program is to perform a close inspection of the inner disk morphology 
Methods. Images were collected with NACO, the AO-assisted near-IR instrument on the VLT (ESO) which includes two types of 
coronagraphs: classical Lyot masks and phase masks. In this program we took advantage of both types of coronagraphs in two spectral 
bands, H-band for the Lyot mask and Ks-band for the phase mask. The Lyot mask blocks a large central region around the star (radius 
< 0.35") but allows deep integrations and hence good signal-to-noise ratio at large distances, while the phase mask allows imaging at 
very close separation (down to ~ 0.15" in theory) but conversely is more sensitive to residual aberrations. In addition, we simulated 
an extended object to understand the limitations in deconvolution of coronagraphic images. 

Results. The reduced coronagraphic images allow us to carefully measure the structures of the debris disk and reveal a number of 
asymmetries of which some were not reported before (position, elevation and thickness of the warp). Our analysis also demonstrates 
the advantage of the phase mask coronagraph to explore the very close environment of stars. In this program, the circumstellar material 
is visible as close as 0.7" (13.5 AU) owing to the phase mask while the Lyot mask generates artifacts which hamper the detection of 
the dust at separations closer than 1.2" (23.2 AU). The point source detection limit is compared to recently published observations 
of a planet candidate. Finally, the simulations show that deconvolution of coronagraphic data may indeed produce artificial patterns 
within the image of a disk. 

Key words, stars: p Pic - stars: pre-main sequence - stars: planetary systems: formation - stars: circumstellar matter - methods: 
observational - techniques: high angular resolution 

1 . Introduction 80 AU) showed that the inner disk (r <25 AU) is depleted by a 

„ , s , . , factor of ~ 100 with respect to the regions of peak density (r « 

5 The young «20Myr) and nearby A-type star /J Pictoris has 80 AU) , Based on new datasets , this d letion was confirmed 

X been the subjec t rfmany investigations since the imaging by b n^etaU ifl997h who proposed a two-component model 

H ; I Smith&Terri el il£84j) of the first circumstellar disk around a tQ int t both the visible/near -infrared and mid-infrared data, 

star. Over the last 25 years, the edge-on dusty disk was observed FrQm moddi | B urrows et al.1 (fl99l derived the presence of 

at many wavelengths, in scattered light and in therma emission, a relativel clear zone within 40 Au In ^ th revealed 

from space and from the ground with the prime goal to under- a Uaj . symmetrical structure , the so-called warp, at a dis- 

stand the planet formation process, and ultimately to find evi- tance <7Q Au wWch should haye disappeared in about l Myr . 

dence tor planetary objects. They interpreted this feature as the gravitational signature of 

The discoverersJ Smith & Ternle| dl984j) .resolved the disk in a lanet object which would also cause the cleari of the 

the visible at projected separations larger than 5.2 (-100 AU, innef disk ions wkh the improvement of Adaptive 0p tics 

adopting the Hipparcos distance of 9.3pc), and extrapolating m (AQ) on large groun d-based telescopes, the warp was 

the observed surface brightness profile to shorter distances they confirmed b M ouil let1taT] lfl9973) in the near-IR using the 

proposed that the disk should be opti callythick below 15 AU 3 6 m te le S c ope at La Silla ( ESQ). To account for the presence 

About ten years latter, [Kalas_&Jewitt| $99$) observed the disk of ^ |Mouillet et al.1 (fl997bh developed a model which 

down to 2.8 (54 AU) and identified morphological and bright- concluded the nce of a lanet on an inclined orbit (3.50) 

ness asymmetries between the north-east and south-west exten- and located in ±c U2Q Au from the star This was 

sions of the edge-on disk. The WFPC2 instrument onboard the alsQ obseryed Jn detail with the HS T/STIS instrument at visi- 

wa i use ^ t0 observe the disk as close as 1.5 (29 AU) We wavelengths dHeap et al.ll2000h . The small Lyot mask and 

and confirms the asymmetrical structure. Mid-infrared imaging ±Q u Qf ±e psp subtraction s allow th e detection 

d Lagage& Pan tin 1994) of the innermost regions of the disk (0- - th . ... as r . . t , star as n 7V , Arf , nrHitia tn | Hp;m Pt a , I 



of the disk as close to the star as 0.75". According to lHeap et aJJ 
Send offprint requests to: A. Boccaletti the Perturbing planet may have a mass of a few times that 
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Angular distance in arcsec Angular distance in arcsec 

Fig. 1. Target to reference flux ratios as a function of the angular separation to the star for the Lyot image (left panel) and the FQPM 
image (right panel). The red and blue curves are the average values measured in angular sectors of 30° and orthogonally to the disk 
midplane (north-west and south-east directions). The vertical bars correspond to the l<x dispersion in these sectors. The adopted 
scaling factors are shown as horizontal dashed lines. 



of J upiter if located clos er than 20 AU. The most recent analy- 
sis dFreistetter et al.ll2007l) based on numerical simulations con- 
strained by the presence of planetesimal belts inferred from the 
observed temporal variations of metallic lines in the j3 Pictoris 
spectrum suggests a planet of 2 My at 12 AU. Note that high res- 
olution spectra were also used to searc h for radial velocity vari- 
ations of ft Pictoris( Gall and et al.ll2QQ6b . The available data al- 
lowed them to detect pulsations (also seen in photometry) and to 
constrain the presence of planets with small periods (e.g. a few 
days). Longer periods planets are still poorly constrained with 
radial velocity data, with for example, at 1 AU, a limit of 9 My. 
While this paper was being review ed, a candidate planet was dis- 
covered by lLagrange et al.1 (120081) taking advantage of a reduced 
contrast in L' band images. With an inferred mass of 8 My, at a 
proj ected distance of 8 AU it is significantly more massive than 
the iFreistetter et al.l (120071) predictions although more measure- 
ments would be needed to definitely rule on the mass of this 
candidate planet and confirm the companionship. 

A detailled d ynamical model of the J3 Pictoris dust disk 
was developed by lAugereau et al.l (120011) who simulate a disk 
of planetesimals and dust particles pertubed by a giant planet 
on an inclined orbit with respect to the disk midplane. The 
Kal as & Jewittl (1 19951) large scale asymmetries (at hundreds of 
AUs) and the warp observed in scattered light are well repro- 
duced assuming a planet located at 10 AU about 10 3 times less 
massive than the central star. The inferred surface density for the 
disk corresponds to a broad annulus peaked around 80-1 10 AU 
and smoothly declining inward within about 70 AU. Although 
the far-infrared thermal emission is correcty predicted by this 
model, it is shown that most of the mid-infrared emission is very 
likely produced by an additional hot population of small grains 
close to th e star and not model ed. Recent high-resolution mid- 
IR images ( Tel esco et al.ll2005l) of the inner disk obtained with 
Gemini have revealed the presence of an asymmetric structure 
located at about 50 AU and which may correspond to a clump 
of particles differing in temperature and/or in size to the other 
particles in the disk. This observation supports the commonly 
shared hypothesis that the disk is replenished by short-lived dust 
originating from collisions of planetesimals but the mechanisms 
that can produce a clump of that size are still unclear. 



The latest ima ges of the (3 Pictoris disk were obtained by 
Golimowski et al. ( 2006) in the visible with ACS on the HST. 
Despite the very high sensitivity of the ACS, the large opaque 
Lyot mask does not allow them to probe the disk morphology 
at an angular separation less than about 1.5" to the star. Using 
a specific data reduction process to emphasize high spatial fre- 
quency structures in the disk, they suggest that the warp is in 
fact the result of a blend of a main disk with a much fainter, in- 
clined secondary disk which can be separated from the main one 
from 80 AU up to about 150 AU in ACS images. However, the 
method used to separate the two disks is arguable (see Sec. 14. II 
for a more detailled discussion) and it is therefore not clear if this 
secondary disk is a fully distinct component or or iginates in the 
main, twisted disk as suppor t ed by t he work of [Mouillet et al. 
(Il997bl) and lAugereau et all (l200ll) . iGolimowski etaD ([2006) 
also showed an interesting color dependence of the disk in scat- 
tered light from which they inferred a typical minimum grain 
size of about 0.2 jim where the planetesimals lie (r < 120 AU). 



Improving the angular resolution in the first few tens of AU 
from ft Pictoris is now clearly mandatory to access the planet 
forming regions. Optimal conditions are met either with large 
ground-based telescopes (providing the wavefront is corrected 
with high-order AO systems) or from space with the HST. In 
this paper, we present new coronagraphic images obtained with 
an 8 m telescope at the European Southern Observatory in the 
near-IR. Although the sensitivity is not sufficient to detect gi- 
ant planets, the coronagraph we used is potentially sensitive to 
smaller separations than any previous observations. In section [2 
we describe the observations we carried out and the data reduc- 
tion process. The analysis of the disk in terms of morphology, 
surface brightness profile, vertical flux and position of the mid- 
plane is presented in section We then discuss the hypothesis 
of a secondary disk (section [47Tb and the contamination by PSF 
structures in coronagraphic images (section l4~2b . 



A. Boccaletti et al.: Fine structures in the disk of J3 Pictoris 



3 





Fig. 2. Images of J3 Pictoris observed with the Lyot coronagraphic mask in the H band. A raw image of the star with the 0.7" Lyot 
mask is displayed to the left while the picture to the right reveals the disk after subtraction of the reference star. The field of view 
is 13.3". The diffraction by the spiders are localized in 3 spikes because the sky is calibrated twice during the observation while the 
pupil rotates. 
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Fig. 3. Images of J3 Pictoris observed with the FQPM in the K s filter (left) and after subtraction of the reference star (right). The field 
of view is 13.3". Here, the sky is calibrated once at the end of the observation. 



2. Observations and data reduction 

2.1. Observing strategy 

Coronagraphic observations of j3 Pictoris were obtai ned on Nov. 
5th, 2004 at ESO with NACO dRousset et alJl2003b the IR AO 
system of the VLT. p Pictoris (V=3.86, H=3.54, K=3.53) was 
observed with 2 coronagraphs, a standard 0.7" diameter Lyot 



mask and a four quadrant phase mask (FQPM hereafter), in the 
H and Ks spectral bands, respectively. Seeing conditions and AO 
compensation were good but variable (Strehl ratios -30-50%). 
Although theoretical FWHMs are 42 mas and 56 mas in H and 
Ks filters, we measured respectively 61 mas and 66 mas result- 
ing from both the actual Strehl ratios and the coronagraphic pupil 
stop (10% undersized). 
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Fig. 4. Point source detection limits (expressed in absolute magnitudes) in the H and Ks bands as measured respectively on the 
Lyot (left panel) and FQPM images (right panel) along the disk (angula r sector of +8°). The absolute magnitudes of some planetary 
masses are overplotted for 10 Myr (derived from the DUSTY models of Chabrier et al. 2000). The 8 Mj limit is reached at 0.55" and 
0.35" for the H and Ks bands. The dashed areas have transmission lower than 50% due to the coronagraphs (inner working angle). 
100% transmission is reached at 0.6" and 0.17" respectively for the Lyot and the FQPM. 




Fig. 5. Magnified image of the central 3" area of Fig. [2] (right). 
The lower circl e shows the expected p osition of the candidate 
planet found by lLagrange et al.l (120081) and the upper circle in- 
dicates the suspicious point-like pattern described in sec. 12.31 
Circles are twice the FWHMs. The dotted line shows the disk 
midplane at PA=29.5° (see section[3]). 



Th e principle of the FQPM is discussed in iRouan etaD 
(2000). It uses a Ti phase shift in the focal plane resulting in a 
destructive interference inside the geometric pupil of a point- 
like object centered on the mask. A monochromatic FQPM 
operating at 2.15/mi was implemented in NACO on August 
2003. Performance as sessment of the FQPM is reported in 
iBoccaletti et al.l (120041) . This coronagraph has no obscuration 



in the focal plane and therefore is able to reach the theoretical 
diffraction limit of the telescope. As a drawback, in actual ob- 
servations (Strehl ratio <50%) the attenuation of the PSF peak 
is less than for a Lyot coronagraph (a factor of ~ 10 instead 
of ~ 200 depending on the level of the AO compensation). The 
chromatis m of the FQPM i s not a limitation on NACO as demon- 
strated in IBoccaletti et al] (120041) . Therefore, we decided to take 
advantage of both coronagraphs to get (i) a good signal to noise 
in the disk at large angular separations using the Lyot corona- 
graph and (ii) to characterize the circumstellar dust at very close 
separations with the FQPM. 

The reference star HD 45291 (V=5.97, H=3.84, K=3.77), 
which was observed right after the target star, was chosen to 
have the same parallactic angle as Pictoris in order to preserve 
the same pupil orientation. This way, the diffraction spikes of 
the reference star have about the same orientation as those of 
Pictoris with respect to the field of view. The reference star has 
the same declination but a difference in right ascension which 
depends on the integration time (36 minutes here). We observed 
(3 Pictoris for a total of 800 s in the H band (Lyot coronagraph), 
and 700 s with the FQPM in the Ks band (but with a transmis- 
sion of only 10% to avoid saturation of individual frames), and 
similarly for the reference star. In addition, the field of view was 
rotated by 45° when the FQPM is used to place the disk in a 
different direction than that of the phase transitions of the phase 
mask. The pixel sampling is 13 mas/pixel. 

2.2. Data reduction 

During the observations, the individual exposures (DIT = 1 s) are 
combined in long exposure images of 100 s each (NDIT = 100). 
Cosmetic reduction is applied on each image including flat field- 
ing, bad pixel correction and sky subtraction. Images are then 
re-centered and co-added independently for the target star and 
the reference star. 

Our procedure to remove the stellar residuals (unattenu- 
ated by the coronagraph) using a reference star is similar for 
the Lyot and the FQP M i mages. As described f or example in 
lAugereau etaD (1 1999b and lBoccaletti et al.l (120031) . the image of 
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Fig. 6. Comparison of J3 Pictoris images and contour lines. The contours on the Lyot image (top) correspond to 12, 13, 14 and 15 
mag/arcsec 2 while on the FQPM image (middle) the lines stand for 12 and 12.5 mag/arcsec 2 . These contours are obviously shaped 
by the numerical sector mask especially in the Lyot image (top). The image at the bottom shows again the FQPM image overplotted 
with the contours of the Lyot image. The spine is shown as a dashed line. 



the target is first divided by that of the reference (re-centered at a 
subpixel scale) to provide a map of the intensity ratio. The scal- 
ing factor is measured at a radial distance where the intensity 
ratio map is almost flat (p = 1.5-5" on the Lyot image and 
p = 0.2 - 1.5" on the FQPM image) and is azimuthaly averaged 
in regions free of circumstellar material (perpendicularly to the 
disk midplane, in the north-west and south-east directions). The 
intensity ratio profiles are displayed in Fig.Q] 

The reference star image is then normalized to the scaling 
factor and subtracted from the target star image. After this pro- 
cess, the Lyot subtracted image still has a diffraction pattern 
originating from the high frequency zonal polishing defects of 
the primary mirr or ([2]). Similar pa tterns are observed in HST im- 
ages of AU Mic (iKrist et al.ll2005l) for instance. This feature was 
estimated in several sectors around the star and an azimuthaly 
averaged pattern was subtracted. However, a residual zonal pat- 
tern remains in the subtracted image and produces a modulation 
of the disk intensity in some particular regions. Nevertheless, 
the disk is clearly visible all across the field of view. The raw 
coronagraphic images and the subtracted images are displayed 
in Fig.[2]and in Fig.[3]for the Lyot and the FQPM. 

For the photometric calibration, out-of-mask images of the 
star were collected, but with a different set-up than corona- 



graphic images. Indeed, these calibration data are obtained with 
a neutral density filter (attenuation of the intensity by factors of 
80 in H-band and 90 in Ks-band) and of a full aperture stop (in- 
stead of a 10% undersized unobstructed circular stop). The pro- 
cedu re for photometric ca libration of such data is fully described 
in iBoccaletti et al.l (120081) . A detailed estimation of photomet- 
ric uncertainties accounting for the differences between out-of- 
mask and on-mask data leads to error bars of 0.1 6 mag in both 
the H and Ks filters. The photometric error budget is actually 
driven by the flux extraction method of a point source (0.13 mag) 
rather than the s caling factors. The sou rces of errors being iden- 
tical to those in IBoccaletti et al.l ([2008), we will adopt the same 
photometric uncertainties. 

2.3. Point source detection limit 

Finally, we measured the point- source detection limits on the 
images of Fig.|2]and Fig.[3j The 5-sigma contrast is calculated as 
the standard deviation of all pixels located at a given radius and 
that for each radius. Once converted to absolute magnitude, the 
result displayed in Fig.|4] shows the detection limit for the H and 
Ks bands along the disk. Note that the 5-sigma detection level is 
taken here for convenience but does not rigorously correspond to 
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a confidence level as usually defined for a Gaussian distribution. 
A detailed study of confidence level in th e presence of diffraction 
residuals is given in lMarois et al.l (120081) . 

The FQPM image (Fig. [3]) has a significantly worse detec- 
tion limit in the background dominated region (>2") owing to 
the 10% transmission. However, the complementarity is clearly 
demonstrated in the limit of detection. The expected absolute 
magnitudes of giant planets (according to the DUSTY model of 
IChabrier et ai1l2000l) for 5, 8, 10 and 20 Mj were overplotted in 
Fig.ffl 

At 0.4", the 4QPM has the ability to pick out an object fainter 
than 8Mj while this level of sensitivity is achie ved at 0.5" in 
the Ly ot image. The candidate planet reported by lLagrange et al.l 
(2008) in Nov. 2003 data should have been detected at least in 
the K band image if locate d at the same p o sition . However, our 
observations and those of lLagrange et al.l (120081) are separated 
by nearly 1 year and hence orbital motion is to be expected. 

The projected distance of the planet candidate was 8AU 
(411 ± 8 mas) on Nov. 2003. We therefore re-analyze our data 
to look around this position but found no firm evidence of a 
companion. However, a suspicious point-like pattern is visi- 
ble in the Lyot image at 511 ± 18 mas i.e. 107 ± 18 mas away 
fro m the discovery position (Fig. [5]). It differs significantly from 
the lLagrange et al.l (120081) observations. This pattern intensity is 
compatible with the DUSTY model (m H « 13.5 - 14). But no 
counterpart is detectable in the 4QPM image (mKQimif) ~ 13.4) 
while the DUSTY model does not predict an equivalent flux in 
the H and Ks bands. Given the proper motion of ft Pictoris (about 
82 mas/yr to the north) a backgro und star would appe ar to move 
closer to the star between 2003 (lLagrange et al.ll2008l) and 2004 
(this paper). The source, if real in both epochs of data, moves far- 
ther away from the star between 2003 and 2004. Despite a lack 
of constraints on the orbital parameters, a projected orbital mo- 
tion of 107 mas can be consistent with an actual separation larger 
than 8 AU. However, the position angle of the s uspicious point- 
like p attern differs by 4.5° with respect to the lLagrange et al.l 
(2008) observations which appears inconsistent with an orbit 
aligned with the disk. We therefore favor a false positive detec- 
tion as an interpretation of the presence of this pattern. 

Two hypotheses can be drawn from our 4QPM detection 
limit at Ks: 

1/ The physical separation in 2004 is less than 8 AU and then 
the planet candidate is at an angular separation less than 0.4" 
which may account for the non detection in our data. In that case 
the motion is greater than 50 mas. 

2/ The physical separation in 2004 is greater or equal 
to 8AU and the planet candidate is fainter than the model 
prediction by almost 1 mag at Ks (M K (limit) « 12). Parameters 
of the system (for instance age) must be reviewed and other 
models could be considered. 



It is difficult with the available material to identify the most 
plausible situation and we will have to wait for a definitive con- 
firmation of the companionship and then carry out additional 
observations in H and Ks with a more effective rejection of the 
speckles. Orbital motion, if observable, will be valuable to derive 
the actual mass of the candidate planet and hence will provide a 
calibration of evolutionary models. 



3. Characterization of the dusty disk 

3.1. Morphology of the disk 

The J3 Pictoris disk is well detected in scattered light in the 
Lyot image (Fig.[6j top panel) up to 7.5" on the south-west 
(SW) side and 7" on the north-east (NE) side (actually limited 
by the detector field of view). At such distances, the surface 
brightness is lower than 15/,mag/arcsec 2 (see section [3T2l) . To 
avoid confusion from diffraction residuals either from the tele- 
scope or the coronagraph, most of the radial features visible in 
Fig. [2] have been masked by selecting a sector of about 30° in 
the Lyot image about the position angle (PA) of the disk (PA 
= 29.5° ± 0.5°, measured from the north and counter-clockwise, 
see sec tion 13. 3b. The brightne ss asymmetry observed in the vis- 
ible by Kalas & Jewittl (119951) is also seen in the near-IR (The 
NE side is brighter than the SW side). The warp is also seen 
when comparing the contour lines to the position of the mid- 
plane (dashed line in Fig. [6]). The largest deviation from the mid- 
plane is observed in between 3" and 4" corresponding to a physi- 
cal distanc e of 59-78 AU, in agreement with the measurement of 
HST/STIS dHeap et al.l2000|) . The innermost regions do not look 
aligned on the midplane but these regions are highly dominated 
by residual diffraction patterns, as explained in the next section. 
Indeed, unlike HST/STIS images, the image of the disk is not 
smooth but is perturbed by the radial diffraction pattern from the 
primary mirror. This diffraction residual produces ringed struc- 
tures that we do not consider as real. 

The FQPM image is displayed in the middle panel of Fig.[6| 
(a sector of 45° about the midplane was selected to avoid confu- 
sion with instrumental radial features). To better distinguish the 
large scale structures from the background noise we applied a 
Gaussian filtering to the image shown in Fig. [3] (right). As ex- 
plained in section [2 the signal to noise is worse than in the 
Lyot image, despite similar integration times, since the trans- 
mission is 10 times lower. The disk is detected down to a level 
of 12.5 mag/arcsec 2 at a radius of about 2". Although the central 
field is not obscured as in the Lyot image, a discontinuity at 0.7" 
on both sides suggests that the patterns located below 0.7" are 
actually diffraction residuals. This is confirmed by the measure- 
ment of the surface brightness in section [3T2l The warp is obvi- 
ously not detected due to a lack of sensitivity at distances larger 
than 2". However, a very pronounced pattern resembling a knot 
at about 1" is clearly visible on the SW side. It was not seen in 
the Lyot image owing to the presence of a diffraction spike and 
its origin remains unexplaine d. This featur e surprisingly resem- 
bles the clump observed by iTelesco et alJ (120051) in the mid-IR 
although not located at the same position (52 AU in the SW side), 
while we are probing closer separations (1" is equiv alent to about 
20 AU ). Ot her clumps (beyond 50 A U) reported bv lWahhai et al.l 
( 2003|) and lWeinberger et a"D (120031) in the mid-IR are interpreted 
as the projection of rings since they are symmetrical. Such rings 
are also observ ed at much l arger separations (>500AU) in the 
visible by lKalas et al.l (l2000l) . However, in our near-IR data, the 
clump has no symmetrical couterpart in the NE side and the ori- 
gin must be found elsewhere than in the presence of rings. The 
refinement of the planet candidate orbit (lLagrange et al] 120081) 
will certainly allows to perform more accurate dynamical sim- 
ulations and possibly be helpful for the interpretation of such 
patterns. 

A smooth contour of the Lyot image is overplotted on the 
FQPM image. The overlap is only partial because the usable field 
is different for each coronagraph (positions of the spikes and 
coronagraph signatures are different). 
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Fig. 7. Surface brightness of the disk (in mag/arcsec 2 ) versus the radial separation for the Lyot (left) and for the FQPM (right). The 
integration times are respectively 800 s in the H band for the Lyot and 700 s in the Ks band with a transmission of 10% for the 
FQPM. The solid lines show the radial profile of the raw image (before subtraction of a reference). The dashed and dash-dotted 
lines stand for the NE and SW sides, averaged over a 16° angular sector. The dotted lines give the noise level as measured at 90° 
from the disk axis. The intensity ranges (1-sigma) in these sectors are displayed as grey shades. 
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Fig. 8. Left: Departure of the disk spine with respect to the midplane as measured on the smoothed Lyot image. Right: Thickness 
of the disk along the midplane. The FWHM of the two sides are overplotted to show the asymmetry. For these 2 plots, parabolic 
fits were used to derive positions of the asymmetries. Error bars are obtained from the same data with several smoothing factors 
relevant to the measurement of such spatial scales. 



3.2. Surface brightness 

The surface brightness (SB) of the disk was measured along the 
NE (PA = 29.5°) and SW (PA = 209.5°) sides and averaged 
over a ±8° sector (avoiding too many diffraction residuals). The 
SB radial profiles derived from the Lyot and the FQPM images 
are shown in Fig. [7] (left and right panels, respectively). The 1 cr 
intensity variation versus radius is displayed with grey shades. 
The surface brightness of the disk is compared to the brightness 
measured in a similar sector but perpendicular to the disk axis 
(a region free of circumstellar material) to estimate an equiva- 
lent noise level. The modulus of this level is plotted in Fig. [7] as 
a dotted line. In the Lyot image, the region closer than -1.2" 
(23.2 AU) is totally dominated by diffraction residuals although 
the coronagraphic mask is only 0.35" in radius. In the FQPM 
image, the field is not obscured in the center but the diffraction 



residuals are clearly stronger than the disk brightness for sepa- 
rations closer than -0.7" (13.5 AU). At large distances (>2.2"), 
the FQPM image is rapidly dominated by the background noise, 
as expected. 

On the NE side, the two profiles are similar within the inten- 
sity range (grey shades), while on the SW side a small difference 
appears near 1.5". However, it is difficult to interpret this differ- 
ence as a color effect because the coronagraphs have different 
signatures and the overlap of the detected regions is only partial. 
Also, it is not sufficient t o constrain the grain siz e distribution as 
we did for HD 141569 in lBoccaletti et al.1 d2QQ3h . 

In both the Lyot and FQPM cases, the diffraction residuals in 
the center can probably be improved with a more careful match- 
ing of the star and reference images during the observation (bet- 
ter centering and better correction of the low order aberrations 
are needed). The advantage of the FQPM over the Lyot coro- 
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Fig. 9. Vertically integrated flux as a function of the angular separation to the star for the Lyot (left) and the FQPM (right). Dotted 
lines correspond to measurements at 90° from the disk axis. 



nagraph, although less sensitive in the background dominated 
region, is clearly seen here and allows us to detect the disk at 
a separation which is almost half the size than on the Lyot im- 
age. This compares to the best distance reported in the literature 
and achieved with the H ST/STIS instrume nt, where the disk is 
detected down to 0.75" ( iHeap et aLll2000|) . However, the exact 
signal-to-noise ratio at t hese distances can not be estimated on 
the STIS images because IHeap et al.l (l200Qh do not compare the 
average brightness of the disk with the noise level. 

3.3. Spine and thickness 

The position angle of disk midplane is measured on the Lyot 
image by a line passing through the star and going from the edge 
of the NE side of the disk to that of the SW side. The disk PA is 
found to be 29.5° ± 0.5°, consistent wit h previous measurem ents 
(between 30.1° and 31.4° according to iKalas & Jewittll 19951 for 
separations larger than ~ 3"). 

In the inner disk regions, the known warp makes the disk 
not perfectly aligned with the midp lane, but instead it ap pears 
twisted as shown for instance by Mo uillet et al.l (fl997b) and 
Hea p et all (|2000). Along cuts orthogonal to the midplane, the 
distance of the peak brightness position to the midplane defines 
the disk spine which is displayed in the left panel of Fig. [8] for 
our Lyot image and for a separation larger than 1.2". This figure 
was obtained in two steps. First, the image of the disk was con- 
volved with a Gaussian function to account for the sole large- 
scale structures. Then, we measured the altitude of the peak 
brightness with respect to the midplane as a function of the dis- 
tance from the star using two different methods: one consists 
of directly measuring the distance of the maximum intensity to 
the disk midplane, while in the second case, the spine position 
is obtained from a Gaussian fitting to the vertical profile. Both 
methods yield very similar results. 

According to Fig. [8] (left panel), the position of the warp in 
the SW and NE sides appears to be not symmetrical with respect 
to the star. The warp extends to 3" on both sides but peaks at 
about 3.3" (65 AU) on the SW side and 4.0" (78 AU) on the NE 
side as measured with parabolic fits (overplotted on Fig. [8]). This 
shift cannot be attributed to a mis-alignment of the star behind 
the mask, since the centering accuracy is better than 0.1". The 
difference in the position of the warp with respect to the star is 



not reported in previous observations. It is reminiscent o f the 
disk offset discovered in Fomalhaut by Kalas et al. (2005) and 
to a lesser extent in the disk around HD 141569 ( Boccaletti et al.l 
2003; Mouillet et al. 2001). The amplitude of this offset is phys- 
ically almost identical in the cases of Fomalhaut and J3 Pictoris 
(about 15 AU), although the ft Pictoris offset might be larger due 
to projection effects. It has been suggested that a planetary com- 
panion on an eccentric orbit (e = 02) co uld create such asym- 
metries (e.g. HR4796. IWvatt et al.|[l999l) . However, it is noted 
that the mass of the pertubating planet is not very well con- 
strained since several combinations of mass and semi-major axis 
are able to reproduce the o bservations. The r ecent detection of 
a planet around Fomalhaut (IKalas et al.ll2008l) confirms the rela- 
tionship between the perturber and the eccentricity of the ring. 
Th e most stringent study in the case of ft Pictoris is presented 
by iFreistetter et al] (120071) . To reproduce the planetesimal belt 
positions inferred from high-resolution mid-IR spectroscopy, a 
system with three planets is necessary, one of them being signif- 
icantly more massive (2^ 5 Mj) and located at 12 AU A rev ised 
value of the J3 Pictoris age ( 12 Myr, |Z uckerma n et al]|200lh to- 
gether with the reasoning of IHeap et al] (l20QQh and an averaged 
position of 7 1 AU of the warp yield a mass of about 2 Mj, sim- 
ilar to that of IFreistetter et al] (120071) . Our detection limit pre- 
sented in Figgis able to place an upper mass limit of 5-6 Mj at 
10 AU for an hy pothetic planet of 10 M yr. This is identical to the 
upper limit of IFreistetter et al] (120071) and therefore not relevant 
to improve the mass c onstraint while the planet candidate from 
lLagrange et al] (120081) is actually more massive (8 My). 

The altitude of the maximum deviation with respect to the 
midplane reaches 0.10" (1 96 AU) in t he SW and 0.07" (1.37 
AU) in the NE extension. IHeap et all (l200Qh also observed a 
higher altitude in the SW side although they did not discuss 
this point. Another deviation from the midplane is also visible 
at closer separations and is located at 1.1" and 1.5" respectively 
on the SW and NE sides. This inverse warp is likely an artifact 
since we do not detect it on the FQPM image which is more ac- 
curate than the Lyot image in this region. We suspect that the 
confusion near the center between circumstellar material and 
diffraction residuals from the on-axis star is responsible for such 
patterns which are enhanced by the use of a Lyot coronagraph . 
It might be the same kind of artifact that lead lHeap et al] (2000) 
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and lGolimowski et all (120061) to conclude that the two extensions 
in the inner disk do not intersect at the star. 

The Gaussian fit also yields the full width at half maximum 
(FWHM) of the disk along the midplane. The maximum thick- 
ness occurs at about 4.5" and 4.0" from the center respectively 
for the SW and NE sides (Fig. [8] right). This asymmetry is op- 
posite to that of the warp position. In addition, the SW side is 
marginally thinner than the NE side (at 1- sigma) a character - 
istic also reported at mid IR wavelengths (iPantin et al.lll997l) . 
The avera ge thicknes s is ab out 0.9" (18 A U), in perfec t agree- 
ment with lHeap etaD (l20QQh and again with lPantin et al.1 (1 19971) . 
There is a consensus fo r this parameter at visible , near IR and 
mid IR wavelengths. In iGolimowski et all (|2006), the SW side 
appears thinner although error bars are missing to confirm this 
point. However, they claim that actual thicknesses are 30% less 
after PSF deconvolution. However, the narrow HST/ACS PSF 
(FWHM=0.05") cannot increase the apparent size of a pattern 
which is more than 15 times larger. With a simple numerical sim- 
ulation, we measured that a disk of which the FWHM is 15 times 
larger than the PSF size undergoes an increase of only 0.5% in 
size once convolved with the PSF. Clearly, the impact of PSF- 
deconvolution in coronagraphic imaging is a very critical issue 
and will be discussed in section |4~2l 

A summary of the asymmetries we measured inside the disk 
is reported in Table [TJ 

3.4. Vertically integrated flux and normal optical thickness 

We finally estimate the vertically integrated flux which is mea- 
sured across the disk thickness (z-axis) and shown as a function 
of the radial distance to the star in Fig. [9] This flux is obtained 
by integrating the disk surface brightness over a Az = 1" region 
centered on the disk midplane, and normalized to the integrating 
area. Residual patterns about the center are avoided with an ap- 
propriate mask which is adapted to each coronagraph. The ver- 
tically integrated flux is almost symmetrical for the NE and SW 
sides in the Lyot image. This is not the case with the FQPM. 
Again, the advantage of the Lyot is clearly seen at large separa- 
tions (>2"), but the FQPM image is still useful to rule out diffrac- 
tion artifacts at separations closer than 2" and down to 0.7". The 
same measurement performed at 90° from edge-on indicates the 
regions where we can be confident in the value of the vertically 
integrated flux. As we took the modulus of the residuals (to avoid 
negative values), a numerical effect makes the equivalent noise 
level sometimes larger than the flux measured inside the disk. 
This is particularly pronounced with the FQPM. 

The disk being seen edge-on, the observed surface bright- 
ness profile at any given projected distance r pro j is the sum of 
the brightness contributions, modulated by the scattering phase 
function, of all the grains physically located at distances r > r pro j 
along the considered line of sight. Deriving the shape of the sur- 
face density profile from an edge-on disk brightness profile re- 
qu ires the use of a deprojecti on method, such the one employed 
by lAugereau & Beustl (|2006) for AU Mic or the numerical one 
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Fig. 10. cr^ ca Z(r) dust disk profiles inferred from the inversion 
of the Lyot vertically integrated SW and NE profiles diplayed 
in Fig. [9l and assuming 3 different scattering asymmetry factors: 
\g\ = (solid line), \g\ = 0.25 (dotted line) and \g\ = 0.5 (dashed 
line). The gray solid line shows the planet esimal pr ofile (normal- 
ized to 0.0045 at its peak position) inferred by lAugereau et al.l 
(2001) for comparison. 



described in IPantin et aD (1 19971) for p Pictorisin this approach, 
the product of the dust surface density profile E(r) and the mean 
scattering cross section of the grains <x sca is reconstructed from 
the outer disk edge down to the closest distance. For the in- 
version process to work properly, smooth profiles that extend 
to large distances where the disk brightness vanishes are re- 
quired. It is therefore not well adapted to the FQPM profiles, 
but more suited to the Lyot surface brightness profiles which 
cover a much larger range of distances. Before performing the 
inversion, the Lyot profiles have nevertheless been extrapolated 
beyond 140 AU with a r -4 5 ra dial profile and a Gauss ian vertical 
profile (FWHM(r) taken from lGolimowski et aT 2006), and then 
smoothed over 10 pixels. 

Examples of depr ojected normal disk pro files using the in- 
version technique of lAugereau & Beustl (120061) are displayed in 
Fig-Hoi for the SW and NE sides. The figure shows the prod- 
uct of the mean scattering cross section of the grains in the 
H-band <x^ a by the dust disk surface density Z(r) for different 
scattering asymmetry factors \g\ between (isotropic scattering) 
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Fig. 11. Figures identical to Fig. [6] (top and middle) but after deconvolution and division by the original, non-deconvolved image. 
Two dashed lines at 5° and -7° from the main disk trace the observed radial structures in the Lyot image. 



and 0.5. The <x^ a X(r) profiles are equivalent to normal dust disk 
optical thickness profiles multiplied by the mean albedo of the 
grains in the H-band. The J3 Pictoris disk appears 2 to 4 times less 
dense at the peak position than the 2.5 times more distant annu- 
lus about the F5/F6 V star HD 181327 that peaks around 87 AU 
(ISchneider et alJl2006h . 

Beyond ~ 150 AU, the calculated profiles follow a r" 2 5 
power law as expected from the extrapolated surface bright- 
ness profiles at these distances. The SW cr^ ca L(r) profiles show 
a peak, whose exact position varies with \g\, but is situated in 
the range 80-90 AU. The global ring-shape of the SW dust disk 
is consistent with the normalized (to allow comparison ) plan- 
etesimal surface density inferred by lAugereau et al.1 ([2001) (gray 
thick line in Fig.fTOb. The dust profiles in Fig. [Kj appear more ex- 
tended than the theoretical planetesimal belt because the radia- 
tion pressure force causes the grains produced th rough collisions 
to po pulate the outer disk regions as shown in lAugereau et al.l 
(2001). The NE cr£ a E(r) profiles below 100 AU do not mirror the 
SW profiles suggesting a non-axisymmetrical structure. The NE 
profiles indeed appear flatter than the SW profiles, and the peak 
at around 90 AU only becomes visible for very anisotropic scat- 
terers (large \g\ values). It is also noticeable that the observations 
in the mid-infrared deduced surface density profiles (iPantin et al.l 
1 19971) are consistent with |g|=0.5 profiles estimated here. 

Although the exact shape of the disk cannot be derived 
uniquely (because of uncertainties in the observations and in 
the scattering properties in particular), both the NE and SW 
cr^ ca L(r) profiles point toward the presence of dust material in- 
side the main 80-100 AU planetesimal belt , as previously sug- 
gested for example by mid-IR observations (iTelesco et al.H2005l 
and ref. therein), and by the present FQPM images. Overall, the 
results tend to support the idea of a complex, structured and 
asymmetric system ins ide the apparent warp location, which led 
iFreistetter et a l. ( 2007) to conclude the presence of several plan- 
etary companions. 



4. On the warp 

The warp of J3 Pictoris has been the subject of many inves- 
tigations and hypotheses esp ecially regarding the presence of 
pertubing bodies like planets (iMouillet et al.lll997bl; iHeap et al 
200(j lAugereau et al.ll200lt IFre istetter et al. 2007). Heap et al. 
(l2000h and later iGolimowski et al.l (120061) have interpreted the 
warp as the blend of two distinct disks at slightly different posi- 
tion angles instead of a single d istorted structure. To em phasize 
the presumed secondary disk, IGolimowski et al.l (120061) used a 
reduction process that we reproduce in Sec. 14. II and critically 
discuss in Sec. 14. 21 

4.1. The secondary disk 

The coronagraphic images presented in iGolimowski et al.l 
( 2006) were obtained with the HST/ACS instrument and classi- 
cally subtracted with the image of a reference star. In this paper 
we adopted a similar approach for the data reduction (Sec. l2.2l) . 
This processing is extremely efficient in the case of the HST 
owing to the stability of the PSF compared to that of ground- 
based telescopes. Then, the HST/ACS subtracted coronagraphic 
image was deconvolved with a st andard Lucy Richard son algo- 
rithm taking a synthetic Tiny Tim dKrist & Hook 2004) image as 
a PSF. 

In classical imaging, the image results from the convolution 
of the object by the PSF. Deconvolution assumes that the PSF 
is stationary (independent of its position in the field). This is 
no longer the case in coronagraphy since the i mage variatio n 
is strongly non-linear especially near the center (lMalbetlll996l) . 
The term PSF is no longer applicable. In a coronagraph, the con- 
volution occurs in the pupil plane downstream of the mask where 
the pupil is c onvolved by the Fourier transform of the mask 
(Malbet 1996). A PSF defined for each point in the image might 
be nec essary in this case for deconvolution. iGolimowski et al.l 
( 2006) have worked around the problem by avoiding the inner 
region inside a 1.5" radius, but the effect of deconvolution on 
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coronagraphic image is not further analyzed. At this stage, the 
warp is still visible as a structur e in the main disk . To im prove 
the visibility of fine structures, iGolimowski et al.l (120061) have 
divided the deconvolved with the non-deconvolved subtracted 
coronagraphic image. This procedure reinforces the high spatial 
frequencies of the images (as in unsharp masking) and reveals 
the warp as two separate disks. This hypothetic secondary disk 
can be distinguished from the main disk beyond 80 AU. Below 
that distance, the intrisic disk thickness prevents any partition of 
the two disks. 

However, this structure is difficult to interpret. It could be 
a secondary disk fully separated from the main disk and cen- 
tered on the star, or could instead originate from the main disk 
as for instance suggested by the inward extrapolation of NE 
second ary disk which intercept s the main NE disk at around 
30 AU (iGolimowski et al.l [20061) . The surface brightness profile 
of the secondary disk as measured by IGolimowski et al.l ( 2006) 
is very reminiscent of the brightness profile of the main disk be- 
yond 125 AU. The latter distance is supposed to represent the 
outer edge of the main planetesimal belt from which the dust 
grains originate that pop ulate the outer disk reg ions due to radia- 
tion pressure forces (e.g. lAugereau et al1l200ll and ref. therein). 
Similarly, the secondary disk could represent the dusty tail of 
an inclined population of planetesimals located within 80 AU. 
It is noteworthy that the extrapolation inward of the secondary 
disk surface brightness measured by Goli mowski et al.l (120061) 
would make it brighter than the main disk below 70 AU, which 
would make the main disk brightness profile much steeper than 
actually observed. The rather flat surface brightness profile ob- 
served below 4" suggests that the secondary disk brightness pro- 
file breaks at around 70-80 AU and becomes flatter inside that 
distance. The reasoning developed for the main disk may thus 
also be valid for the secondary disk, which would then originate 
from parent bodies arranged in the form of an inclined belt with 
an outer edge close to 70-80 AU. It is remarkable that this dis- 
tance matches the maximum extensi on for the warp produce d by 
the inclined planet in the models of Mouill et et al.l ([1997a) and 
lAugereau etaD (120011) . 

Our reduced coronag raphic images were deconvolved as in 
IGolimowski et al.l (120061) except that we were using an actual 
PSF image of J3 Pictoris (out of the coronagraph) and a maxi- 
mum likelihood algorithm for the deconvolution. As a result of 
the deconvolution process, we marginally detect a radial struc- 
ture with an inclination of 5° with respect to the disk midplane 
and consistent with the secondary disk, but also one almost sym- 
metrical structure about the disk midplane at -7° from the main 
disk (Fig.Qj] top panel). It is however difficult to measure the ra- 
dial distance where these structures originate, although they are 
consistent with purely radial patterns intersecting the star unlike 
the secondary disk (at least on the NE side) observed with the 
ACS instrument. Here, we interpret the radial patterns revealed 
by the deconvolution process as PSF radial patterns instead of 
real pattern s. These patterns are th erefore different to those ob- 
served by IGolimowski et all (120061) . Our emphasized FQPM im- 
age is presented in the bottom panel of Fig.[TT]but does not show 
any features. 

4.2. Deconvolution of coronagraphic images: simulation 

To evaluate the significance and the nature of the PSF patterns 
that may appear after deconvolution of a coronagraphic image of 
a circumstellar disk, we produce synthetic HST/ACS-like coro- 
nagraphic observations of an edge-on disk that we afterwards 
deconvolve. Instead of using a Tiny Tim PSF, we produce a syn- 




Fig. 12. Simulated images of (from top to bottom), the disk 
model, the disk convolved with the PSF, the coronagraphic im- 
age subtracted with that of a reference star and the deconvolved 
image. The FOV is 4.5" x 1.3". 



thetic HST PSF image using our own diffraction code that we 
have developed to model the behavior and performance of coro- 
nagraphs, as well as to gauge the potential of several plane t find- 
ing projects (lBoccalettill2004l : lBoccaletti et al .1120051 120061) . 

We consider a 2.4 m diameter telescope observing at 
0.606 yum, with a bandwidth of 33%, and providing an angular 
resolution of 0.052" (similar to the ACS parameters). We gen- 
erate the image of an edge-on circumstellar disk having a 4.5" 
diameter extension and a midplane surface brightness propor- 
tional to r~ L5 , r being the angular distance to the star. The disk 
extends vertically over 0.52" (10 times the angular resolution) 
with a Gaussian profile of FWHM = 0.17" (Fig.El top panel). 
The disk model is absolutely free of structures. Note that some 
of these parameters are arbitrary and are not intended to exactly 
model the J3 Pictoris disk. This object contains 8542 points for 
which we have to calculate the image through a coronagraph (a 
Lyot mask of 1.8" diameter) and for several wavelengths so as to 
produce a realistic coronagraphic polychromatic image resem- 
bling that of the ACS instrument. To simulate realistic images, 
the model accounts for the phase aberrations u sing the measured 
phase map of the HST dKrist & Burrowslfl995l) . and we assume 
similar static aberrations on the reference star but with an addi- 
tional aberration in the form of a defocus (4 nm rms) to account 
for the so-called "PSF breathing" (the HST PSF is slightly evolv- 
ing over time). This amount of defocus is taken arbitrarily as we 
are looking for qualitative results here. Another way to produce 
differences in the target star and reference star images would be 
to apply a different pointing for the two images with respect to 
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Fig. 13. Comparison of the radial intensity profile (left) and the vertical profile of the disk obtained at 1.5" (right) for the several 
steps of the simulation. 



the Lyot mask. The star to disk flux ratio is also arbitrary and is 
selected to have the residuals on the subtracted image of the or- 
der of the disk intensity. We suspect that this is precisely in this 
regime of contrast that deconvolution may generate artifacts. The 
synthetic images is then deconvolved as in Sec. 14. II 

Several images to illustrate the process of image forma- 
tion are displayed in Fig. [12] The two top panels show the disk 
model alone (no star) and the disk convolved with the PSF. The 
two lower panels are coronagraphic images (reference star sub- 
tracted) before and after deconvolution. Several patterns like 
speckles, rings, and radial spikes which were not included in the 
disk model are seen and are clearly enhanced in the deconvolu- 
tion process. For instance, the radial spikes mimic the splitting 
of the disk midplane into two parts but are obviously not real. 
These spikes are aligned here with the disk to depict a worse 
case but could be differently oriented. 

We used these simulated data to measure the averaged radial 
profile and the disk thickness (Fig. [13]), before and after decon- 
volution of the image. Despite the intensity fluctuations resulting 
from the PSF patterns the averaged radial profile is consistent 
with the model except near the mask edges. It is interesting to 
see how the deconvolution enhances the variations of the inten- 
sity profile (Fig. [13] left) especially at distances closer than 1.0- 
1.2". As a confirmatio n of the discussion in Sec.J3.3l and contrary 
to what is proposed in lGolimowski et aD ([2006), the disk thick- 
ness is not affected by the deconvolution as long as it is thicker 
than the angular resolution (Fig. [13] right). In the simulation, the 
FWHM is only 0.17" which is about 3 times the angular reso- 
lution while in the actual data of ft Pictoris, the disk thickness 
is about 0.87" (Fig. [8]) so the PSF blurring should be even less 
visible. The vertical cut shown in Fig. [13] (right) is obtained at 
1.5" (left side of the image). The shape of the blue line (decon- 
volved image) changes significantly at closer separations, so the 
FWHM of the disk must be measured in a region less affected 
by PSF patterns. 

Although the simple simulation presented here reveals the 
production of artifacts, a more general analysis would be neces- 
sary to place quantitative limitations when using deconvolution 
with coronagraphic images of circumstellar disks. We note that 
the simulations presented here d o not rule out the the rel iability 
of the secondary disk shown in iGolimowski et all (120061) . This 
secondary disk is clearly significant with respect to PSF patterns 



arising in deconvolution of coronagraphic images. In particular, 
it is symmetrical and originates in the main disk at a distance 
where PSF patterns are much less preponderant (see Sec 14.11) . 
The PSF patterns that may be confused with disks structures like 
spikes mostly originate from the star. 

5. Conclusion 

This paper reports on the observation of the J3 Pictoris disk in the 
near IR with NACO at the VLT using two coronagraphs in two 
different spectral bands. We have used a careful data reduction 
process to get rid as much as possible of the diffraction residuals 
and artifacts that arise in high contrast imaging with a corona- 
graph on ground-based telescopes. The resulting images reveal 
the dusty disk in preferential regions depending on the type of 
coronagraph: between 1.2" to 7.5" with the Lyot mask and be- 
tween 0.7" to 2.5" for the FQPM image. These two sets of data 
are therefore complementary in terms of spectral bands and an- 
gular separations. 

The Lyot mask coronagraph allows a very good characteri- 
zation of the disk at large distances and especially of the famous 
warp. We have been able to precisely measure the position, el- 
evation and thickness of the warp and we have shown geomet- 
ric and photometric particular asymmetries that could eventually 
been used as inputs to dynamical modeling in order to better con- 
strain the characteristics of the forming planetary system. Some 
of these asymmetries were reported in various spectral bands 
(visible and mid IR). The inversion of the surface brightness pro- 
files have confirmed the presence of an asymmetric main plan- 
etesimal belt peaking at around 80-100 AU responsible for the 
dust material beyond that distance. The inversion further sup- 
ports the presence of significant amounts of additional material 
inside the main planetesimal belt as suggested by other observa- 
tions, depicting a complex planetary system. 

The FQPM image reveals the presence of very close patterns 
very similar to what is observed in the mid-IR, a spectral regime 
that is sensitive to the dust emission. The interpretation of such 
structures would require further analysis possibly involving dy- 
namical models to understand if the origin is similar to those 
observed in the mid-IR. 

To understand the issue of deconvolution in coronagraphic 
imaging we have performed numerical simulations to mimic 
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HST Lyot coronagraphic data. Assuming a disk with no struc- 
ture, this qualitative work shows how the data reduction process 
may emphasize some PSF structures that reinforce the need for 
a thorough analytical study for this type of data. 

These new images of the J3 Pictoris debris disk presented 
here demonstrate the interest of ground-based instruments to 
achieve high angular resolution at longer wavelengths than with 
the HST and motivate for even more accurate observations. The 
sub-arcsec region is now receiving most of the attention with 
the recent discovery of a planet candidate. Constraints on the 
characteristics of this object and possibly others could be in- 
ferred form the spatial distribution of the dust. Therefore, the 
precise registration of the inner structures we have performed if 
complemented with perturbation theory will help to better char- 
acterize thi s planetary sys t em. To confirm the numerical simu- 
lations of iFreistetter et al.l (120071) the study of the disk and its 
structures near 12 AU is essential. This corresponds precisely to 
the limit reached in our data (0.62") and therefore the motivation 
is strong to perform more accurate observations. This is becom- 
ing possible at the VLT. The coronagraphic suite of NACO now 
implements other FQPMs tuned for the H and Ks bands with 
the ability to reach a better attenuation. Additional techniques 
also have been installed (Angular Differential Imaging, Spectral 
Differential Imaging) that may facilitate the calibration of the 
residual starlight and hence the registration of structures even 
closer to the star. 

After more than 20 years of intensive studies, the dusty disk 
around ft Pictoris is still the subject of many questions. Direct 
imaging observations are complementary to spectroscopy and 
photometry and are mandatory to better constrain the spatial dis- 
tribution of the dust, planetesimals and hypothetical planets in 
order to contribute to the understanding of the planetary forma- 
tion phase. 
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